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Abstract. 

We introduce a generalization of the Dobihski relation through which we define a family 
of Bell- type numbers and polynomials. For all these sequences we find the weight 
function of the moment problem and give their generating functions. We provide 
a physical motivation of this extension in the context of the boson normal ordering 
problem and its relation to an extension of the Kerr Hamiltonian. 
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In this note we consider the Dobihski relation and its generalization. This topic 
naturally belongs to the field of combinatorial analysis. The Dobiiiski relation [1] 
was first derived in connection with Bell numbers B{n) = 1,1,2,5,52,203,877,..., 
n = 0,1,2,..., which describe partitions of a set [2], [3]. That remarkable formula 
represents the integer sequence B{n) as an infinite sum of ratios 

oo un 

fc=o ^■ 

Closely related to the Bell numbers are Stirling numbers of the second kind S{n,k), 
k = l...n, and the Bell polynomials defined as 

x) = A;) x^ (2) 

k=l 

related to B{n) by B{n) = B{n,l) = YIk=iS{n,k). For the Bell polynomials the 
Dobihski relation Eq.(l) generalizes to 

oo jLn 

B{n,x)^e-^Y.1J^'- (3) 
fc=o ■ 

These formulas may be derived using either combinatorial or purely analytical methods 
starting from the original interpretation of Bell and Stirling numbers given in 
enumerative combinatorics [4], [5]. Accordingly, Stirling numbers S{n,k) count the 
number of possible partitions of the n element set into k subsets (none of them 
empty) and Bell numbers B{n) count all such partitions. We note that other pictorial 
representations can also be given, e.g. in terms of graphs [6] or rook numbers [7], [8], [9]. 

One may conversely take Eq.(l) (or Eq.(3)) as the definition of the Bell numbers (or 
polynomials). This observation suggests the generalization of these sequences through 
the Dobihski relation. In this note we introduce an extension of Eq.(l) and define the 
family of Bell-type numbers as 

where P{k) and D{k) are any functions of A: = 0, 1,2, ... such that D{k) ^ and the 
above sum converges. Note that conventional Bell numbers are obtained for P{k) = k 
and D{k) = ek\ . 

This generahzation was also pointed out in [10] and [11] in connection with the log- 
normal distribution. Here we focus on the general properties of our proposed definition 
and show that the very specific form of Eq.(4) results in a straightforward solution of 
the moment problem and calculation of the generating functions. We also comment on 
the connection to physics and interpret the sequences so defined in the context of the 
problem of the normal ordering of boson operators. 

Suppose that we want to solve the moment problem [12] for the sequence B{n), i.e. 
we seek a positive weight function W(|/) such that B{n) is its n-th moment 

B{n)^ Jdyy^Wiy). (5) 
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Figure 1. The portion for < y < 5 of the weight function W{y) = X^fe^o fcT 
solving the Stieltjes moment problem for the Bell numbers B{n). Height of the vertical 
lines is proportional to the weight of the Dirac 6 functions. 



At this point wc do not specify the domain of W(|/) or the hmits of the integral. A 
closer look at Eq.(4) yields the following candidate for the weight function 

This is an infinite ensemble of weighted Dirac S functions located at a specific set of 
points {P{k), k = 0,1,2,...} and is called a Dirac comb. If all the weights 1/D{k) 
are positive {D{k) > 0) and normalized to one {J2T=o^/^i^) ~ ^) then Eq.(6) is 
a positive and normalized distribution which is a solution of the moment problem 
of Eq.(5). Whether it corresponds to the Hamburger, Stieltjes or Hausdorff moment 
problem depends on the range of the set {P{k), k = 0,1,2,...}. For example, for 
the sequence of Bell numbers B{n) the weight function W{y) — X^fclo ^^^F^ ^ 
positive and normalized distribution solving the Stieltjes moment problem, see Fig.l. 
A solution of the Hamburger moment problem is generated by the set of restricted Bell 
numbers B-{n) = 1, 0, 1, 1, 4, 11, 41, 162, . . . for n = 0, 1, . . . counting partitions without 
singletons [4]. They satisfy Bj{n) = SfcLo ^i^'' with P(k) = k — 1 and D(k) = e k\; 
the measure is yVj{y) = Y^kLo ^^^^~^r^- the other hand, the well-known Catalan 

numbers C{n) — (^) solutions of the Hausdorff moment problem [14]. 

The specific form of Eq.(4) simplifies the calculation of the generating functions. 
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Taking the exponential generating function, substituting Eq.(4) and changing the 
summation order one obtains 

oo \n oo \P{k) 

Q{\) = T.mzj = Y.-TuTx- (7) 



n=0 



to m ■ 



Evaluation of this series depends on the particular choice of the functions P{k) and 
D{k) and in general the series may be divergent. For the Bell and restricted Bell 
numbers Eq.(7) it can be evaluated easily: G{X) = J2'^=o ^i''^)^ — e"^^"^ and 

Similarly, for the ordinary generating function one gets 

oo oo -1 

g„(A) = EHWA" = E ^(,).(,_p(,),) . (3) 

The same procedure can also be performed for other cases, e.g. for hypergeometric 
generating functions [13]. The choice of the denominator in the generating function 
may depend on P{k), D{k) and the purpose we need it for {e.g. when we need a well 
convergent generating function for analytical calculations). 

In the same manner one could generalize Eq.(3) and define 

= (^) 

The additional variable in the functions P{n,x) and D{n,x) does not pose any 
complication either in the proposed approach to the moment problem or in the 
evaluation of the generating functions. However, wc must observe that in general i3(n, x) 
has an infinite expansion in x and only for particular choices of the functions P{k,x) 
and D{k,x) does it yield polynomials. It is certainly the case for D{k,x) = Me^x"^, 
and P{k, x) — P{k) a polynomial in k. Therefore we define the Bell-type polynomials as 

-B(n,a:) = e-^f;£^x^ (10) 
fc=o ^■ 

As a result the weight function of Eq.(6) takes the form 

>v(x,y) = e-f:^^^^^x^ ai) 

A;=0 

and the exponential generating function of Eq.(7) is 

oo \P{k) 

g{X,x)^e-^^—^x'. (12) 
fc=o 

The case of conventional Bell polynomials is obtained for P{k,x) = k. Consequently 
one gets the positive and normahzed weight function W{y, x) — e"^ J2kLo ^^^r^ 
of the Stjelties moment problem, B{n,x) — dy y''^W{y,x), and the exponential 
generating function G{X,x) = J2^=qB{'>T',x)^ = e^^^^~^^ (see [15]). Analogous 
considerations can be applied to the polynomials generated by Bj[n) leading to 
Bj{n,x) = e--Er=o^^'> and G-{X,x) = En=o B-{n,x)^ = e^(^'-'-'\ 

Introduction of the generalized Bell-type numbers and polynomials through Eqs.(4) 
and (10) is not merely a mathematical definition but it has a firm grounding in physics. 
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We will show that it is related to the solution of the normal ordering problem for a 
general function of the number operator. 

Consider the boson creation operator and annihilation operator a satisfying the 
commutator [a, a^] = 1. Suppose we are given a function of these operators. Its 
normally ordered form is obtained by moving all the creation operators to the left of the 
annihilation operators using the commutation relation. The normal ordering procedure 
is of fundamental importance in quantum mechanical calculations in the coherent state 
representation, the latter defined by the coherent states \z) — e"'^'^/^ tIt l'^)' where 
a^a\n) — n\n), {n\n') — 5n,n' ^-nd a\z) — z\z) [16]. For example, if we take the n-ih. power 
of the number operator a^a the normal ordering procedure gives [17] 

(ata)" = A;) (a^) V. (13) 

It involves Stirling numbers of the second kind S{n, k) and the coherent state matrix 
element yields the Bell polynomial 

{z\{a^aY\z)^B{nM^)- (14) 
Now we consider a general polynomial of the number operator denoted by 

N 

Haia^a) = Yl afe(a^a)' (15) 

k=No 

with some constants a^; A'o and N are the smallest and largest indexes of non- vanishing 
ctfe, respectively. Physically, Tia may be thought of as a generalisation of the Kerr 
Hamiltonian of quantum optics [18]. The n-th power of Tia defines the Stirling-type 
numbers as 

nN 

7f^(ata)]" = ^ 5«(n,fc) (at)'=a^ (16) 
fe=Aro 

and associated Bell-type polynomials (of order nN) are 

nN 

B^in,x)= Sain,k)xK (17) 

k=No 

We will show now that such polynomials defined in the normal ordering problem 
correspond to the Bell-type polynomials introduced in Eq.(lO). 

To this end observe that [a, a''"] = [-D,X] = 1 where D and X are the derivative and 
multiplication operators. We first rewrite Eq.(16) in terms of D and X: 

nN 

[n^{XD)r = E Mn, k) X^D\ (18) 

k=No 

By acting with the r.h.s. of Eq.(18) on one obtains e^Ba{n, x). Action of the l.h.s. on 
is more involved. First we apply it to the monomial x"^ which yields \Ha{XD)Y x"^ = 
(EtjVo«fc^') ^"^ from which [n^{XD)re- = E~=o (E^ato «fe S Allows. 

Combining these two observations we deduce that 

Bo^in, x)=Y: Sain, k) x' = e'^ ^ ^^^^ x\ (19) 

k=No k=0 
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which has the same form as Eq.(lO) for P{k) — Ti.a{k). This gives an interpretation of 
the Bell-type polynomials and numbers in the context of the normal ordering problem. 
The assumption that Tiaix) is a polynomial guarantees that Ba{n, x) is also a polynomial 
in X. Although this additional assumption may be irrelevant in general, as we have 
mentioned above it leads to infinite sequences of Stirling-type numbers. 
Having solved the normal ordering problem and identified the solution in the framework 
of the Dobinski-type relations we may now exploit its advantages. The solution of the 
moment problem of Eqs.(5) and (6) is immediate (depending on the coefficients ak in 
Ha{x)). Also the generating function setting is easily appficable. We just mention 
that the normally ordered exponential of a function of the number operator is the 
exponential generating function of the associated Bell-type polynomials. In the coherent 
state representation it may be written as 

k=Q ^■ 

This and other consequences of the normal ordering of polynomial-type Hamiltonians 
in the number operator will be discussed elsewhere. 

In conclusion we want to emphasize the advantages of introduction of the Bell-type 
numbers and polynomials through generalization of the Dobinski relation of Eqs.(4) 
and (10). It enables a straightforward solution of the moment problem in the form 
of a Dirac comb, see Eqs.(5), (6) and (11). Moreover, calculation of the generating 
functions simpfify considerably in that framework (see Eqs.(7), (8) and (12)). We have 
also pointed out that this generalization has immediate application to the boson normal 
ordering problem. We have interpreted a wide class of Stirling-type numbers as the 
expansion coefficients of normally ordered functions of the number operator. 

Further modifications of the structure of the infinite series of Eqs.(4) and (10) 
may, in general, lead to moment problems with continuous weight functions and will be 
developed elsewhere. 
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